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ABSTRACT 

We present collapse simulations of 100 Mq , turbulent cloud cores threaded by a strong 
magnetic field. During the initial collapse phase filaments are generated which frag- 
ment quickly and form several protostars. Around these protostars Keplerian discs 
with typical sizes of up to 100 AU build up in contrast to previous simulations ne- 
glecting turbulence. We examine three mechanisms potentially responsible for lowering 
the magnetic braking efficiency and therefore allowing for the formation of Keplerian 
discs. Analysing the condensations in which the discs form, we show that the build-up 
of Keplerian discs is neither caused by magnetic flux loss due to turbulent reconnec- 
tion nor by the misalignment of the magnetic field and the angular momentum. It is 
rather a consequence of the turbulent surroundings of the disc which exhibit no coher- 
ent rotation structure while strong local shear flows carry large amounts of angular 
momentum. We suggest that the "magnetic braking catastrophe", i.e. the formation of 
sub-Keplerian discs only, is an artefact of the idealised non-turbulent initial conditions 
and that turbulence provides a natural mechanism to circumvent this problem. 
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1 INTRODUCTION 

In recent years a great number of simulations have been per- 
formed that investigate the formation of protostellar discs 



under the influence of magnetic fields (e.g. 


Allen et al.|2003 


Matsumoto & Tomisaka 


2004; Machida et al.||2005| |Baner- 


jee fc Pudritz|2006| |2007 


|Price & Bate|2007| |Hennebelle & 


Fromang 2008; Hcnncbcl 
2009; |Commergon et al. 


e fc Ciardil:2009 
2010| |Peters et 


Duffin & Pudritz 
al.||2011| ISeifried 


et al.||2011 L In simulations with magnetic field strengths 



rart et al. 2009 Beuther et al. 20101 no rotationally sup 



ported discs were found. As strong magnetic braking is re- 
sponsible for the removal of the angular momentum, this 
problem is also called the "magnetic braking catastrophe". 
The results of these numerical simulations stand in sharp 
contrast to observations, which show that discs are present 



in the earliest stage of protostellar evolution (e.g. J0rgenscn 
|et al.|2009[ ). It is also well known that outflows are the first 
observable signatures of star formation, and that therefore 
discs should already have formed (e.g. Arce et aL]|2007 l 



The inclusion of ambipolar diffusion jMellon fc Li|2009 
Duffi n fc Pudritz|2009 1 also fails to produce Keplerian discs 



in the earliest evolutionary stages. Considering the effect of 
Ohmic dissipation, only very small (~ 10 solar radii) Keple- 
rian discs were found |Dapp & Basu 20fT| |Dapp et al.|2011[), 



unless an unusually high resistivity is used ( Krasnopolsky 
et al.||2010| ). However, recently two mechanisms solving the 



catastrophic magnetic braking problem were proposed, i.e. 



the inclusion of the Hall effect (Krasnopolsky et al. 20111 
and turbulent reconnection ( Santos-Lima et al.|2012 |. 

In this paper we remedy a shortcoming of the simula- 
tions referred to in the beginning namely their lack of tur- 
bulent motions. We present results from a number of sim- 
ulations investigating the possible role of turbulence in re- 
ducing the magnetic braking efficiency and allowing for the 
formation of protostellar discs. 



2 INITIAL CONDITIONS 

We now shortly describe the basic simulation setup. For 
a more detailed description we refer the reader to |Seifried| 
et al. (20111. We simulate the collapse of a 100 Mq molec- 



dseifriedOhs . uni-hamburg.de 



ular cloud core, 0.25 pc in size and embedded in 0.75 pc 
sized cubic simulation box of low-density gas (4.2 x 10 -21 g 
cm -3 ). The density in the core declines as p oc r -1 ' 5 hav- 
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Table 1. Initial conditions of the performed simulations including 
run 2.6-4 without turbulence presented in Seifried ct al,|p011| . 



Run 


/' 


A-ot 


Aurb 


turbulent seed 


P 


2.6-4-A 


2.6 


0.04 


0.04 


A 


5/3 


2.6-4-B 


2.6 


0.04 


0.04 


B 


5/3 


2.6-4-C 


2.6 


0.04 


0.04 


C 


5/3 


2.6-4-poly 


2.6 


0.04 


0.04 


A 


5/3 


2.6-4-A-b 


2.6 


0.04 


0.04 


A 


2 


5. 2-4- A 


5.2 


0.04 


0.04 


A 


5/3 


2.6-4 


2.6 


0.04 










ing a maximum of 2.3 x 10~ 17 g cm" 3 in the centrtQ The 
core is threaded by a magnetic field in the z-direction declin- 
ing radially outwards with R~ ' 75 . It has a strength of 1.3 
mG in the centre corresponding to a mass-to-flux ratio of 
H — 2.6. The core is rotating rigidly around the z-axis with 
a rotational energy normalised to the gravitational energy 
of Piot = 0.04. Additionally, we add a supersonic turbulence 
spectrum with a power-law exponent of p = 5/3. The turbu- 
lent energy is equal to the rotational energy, i.e. /3turb = 0.04, 
corresponding to a turbulent rms-Mach number of ~ 2.5. 

The applied cooling routine ( Banerjee fc Pudritz|2006 1 
takes into account dust cooling, molecular line cooling and 
the effects of optically thick gas. We introduce sink particles 
above a density threshold of p cr it = 1.14- 10 -10 g cm -3 (see 
for details). A maximum grid resolu- 



Federrath et al. 2010 



tion of 1.2 AU is used. The refinement criterion used guar- 
antees that the Jeans length is resolved everywhere with at 
least 8 grid cells although a even higher resolution has been 
suggested (Fe derrath et al.|2011[ ). 

To check whether our results are affected by the ran- 
dom turbulent realisation, we performed two more simula- 
tions with identical initial conditions but different turbulent 
seeds. Furthermore, we performed three more simulations 
with 1) a magnetic field with half the fiducial strength 2) 
a power-law exponent of the turbulence spectrum of p = 2 
and 3) a polytropic cooling to explore the dependency of our 
results on the initial conditions and numerical methods. All 
simulations are listed in Table [l] 



3 RESULTS 

In total we have performed 6 simulations with varying ini- 
tial conditions. In the following we present the results of our 
fiducial run 2.6-4-A in detail. However, we emphasise that 
the results of the remaining runs are qualitatively very sim- 
ilar. Hence, we are confident that the main findings do not 
depend on the randomly chosen, initial turbulence field. 

As we are interested in the properties and the evolu- 
tion of protostellar discs, we restrict our consideration to 
the time after the first sink particle has formed (to). De- 
pending on the simulation, this happens after roughly 15 
- 20 kyr. At this point large filaments have developed in 
which the discs form (for fragmentation and disc formation 
in massive cores without magnetic fields see Banerjee et al.| 
2006 Girichidis et al. 112011 1. From this point on we run the 
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Figure 1. Radial dependence of the rotation (red) and radial 
velocity (green) for the four discs formed first in run 2.6-4-A. 
The black solid line shows the Keplerian velocity v^ ep , the dotted 
lines 50% and 150% of v^ ep . The solid line in the negative velocity 
range shows the free-fall velocity Vg = V2v^ ep . The regions below 
4 AU are affected by resolution effects, therefore they are shaded 
grey to guide the reader's eye. 



simulations for further 10 - 15 kyr. In order to determine 
global disc properties like mass, centre-of-mass (CoM) and 
the angular momentum vector we only consider gas with 
densities larger than 5- 10~ 13 g cm~ 3 around a sink particle. 
From visual inspection of the density isocontours and varia- 
tion of this threshold we found this value to be reasonable. 
Furthermore, it roughly corresponds to the threshold where 
the gas gets optically thick. With respect to the CoM and 
the orientation of the disc we now can calculate the rotation 
velocity and the radial infall velocity u ra d for each grid 
cell less than 20 AU above/below the plane defined by the 
disc. In order to get an impression of the scatter of these 
quantities, we do not azimuthally average the values of v^, 
and Urad- In Fig. [I] we show the radial dependence of v$ 
and D ra d 15 kyr after the formation of the first sink particle 
for the four discs formed first in run 2.6-4-A. Some of the 
discs have already fragmented and contain more than one 
sink. The remaining sinks have either no associated discs 
due to an ejection event caused by many-body interactions 
or have been created only shortly before the end of run 2.6- 
4-A so that their discs are not yet well developed, although 
the velocity profiles reveal Keplerian disc features as well. 
We note that, in general, the angular momentum vectors of 
the discs are well off the z-axis. This demonstrates that the 
discs are created by the local angular momentum associated 
with turbulent motions and not by the overall rotation of 
the cloud core. The orientation of the discs does not vary 
significantly over time as the large-scale structure in which 
they reside does not change much during the time consid- 
ered (~ 10 kyr). To get an impression of whether the discs 
are rotationally supported or not, we also plot the Keplerian 
velocity t>kc P = ( GA ^ star ) 1|/2 in Fig. [lj where G is the grav- 
itational constant and M sta r the mass of all sink particles 
in the disc. As the disc mass (~ 0.1 Mq) is significantly 



1 To avoid unphysically high densities in the interior of the core, 
we cut off the r~ 15 -profile at a radius of 0.0125 pc. 



smaller than the mass of the sinks we easily can neglect it. 
As can be seen, is close to «k cp out to a radius of a few 10 
AU (discs 3 + 4) up to 100 AU (discs 1 + 2) with a scatter 
of about 50% in each direction as indicated by the dotted 
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Figure 3. Mass-to-flux ratio /i (left) and inclination of the mean 
magnetic field to the angular momentum vector of the disc (right) 
in spheres with a radius of 500 AU around the CoM of the discs 
found in run 2.6-4-A. 



Figure 2. Column density in logarithmic scaling for the top-on 
view of disc 1 (top left) and disc 2 (top right) of run 2.6-4-A and 
of the disc in run 2.6-4 without turbulence (bottom). The figures 
are 800 AU in size. 



lines. This is a remarkable result since for previous simula- 
tions of low- and high-mass cores with mass-to-flux ratios 



pi < 10 only sub-Keplerian discs were found (e.g. Allen et al 



2003;|Price & Bate 2007| |Mellon & Li 


2008; |Hennebelle & 


Fromang(2008 


Duffin & Pudritz|2009 


Seifried et al.|2011 


1. 



We emphasise that for the other runs we find qualita- 
tively similar results, i.e. discs with sizes of up to ~ 100 AU 
and masses of the order of 0.1 Mq. The number of discs per 
run varies between 2 and 5. We briefly note that the discs 
presented here drive molecular outflows. Furthermore, in all 
discs v,j, scatters around Ukc P , indicating that this is neither 
a consequence of the specific turbulence seed (runs 2.6-4-B 
and 2.6-4-C) nor of the adopted cooling function (run 2.6-4- 
poly) nor of the power-spectrum exponent (run 2.6-4-b). We 
find that f ra d scatters around and is almost always smaller 
than Vt/y and significantly smaller than the free-fall velocity 
Vg = x/^kcp- This is in strong contrast to the disc in run 
2.6-4 which has the same initial setup as the runs presented 
here except the initial turbulence field ( |Seifried et al.|20"TT I 



The disc was found to be strongly sub-Keplerian with v ra d 
close to vg. The difference becomes particularly clear when 
comparing the top-on view of disc 1 and 2 in run 2.6-4-A 
with that of the disc in run 2.6-4 (Fig. [Sf. 

Why, even in the case of such strongly magnetised cores, 
are Keplerian discs formed? The suppression of Keplerian 
disc formation in previous studies without turbulence is due 



to the very efficient magnetic braking ( Mouschovias & Pa- 



leologou|l980[ ) which removes angular momentum from the 
midplane at a very high rate. Hence, in our runs the mag- 
netic braking efficiency has to be reduced significantly. Two 
possible reasons for this are the loss of magnetic flux in the 
vicinity of the discs or, as proposed recently, a misalign- 
ment of the magnetic field and the angular momentum vec- 



tor ( |Hennebelle fc Ciardi|2009] |Ciardi fc Hennebelle|2010 | 

We first consider the possibility of magnetic flux loss 
in the vicinity of the discs which might be attributed to 
turbulent reconnection (Lazarian & Vishniac 19991. For this 



purpose we calculate the volume-weighted, mean magnetic 
field (B) in a sphere with a radius of r = 500 AU around 
the CoM of each disc. In combination with the sphere mass 
M we obtain the mass-to-flux ratio 



A* = 



M 



0.13 

7u 



(i) 



We plot the time variation of fj, in the left panel of Fig. [3] 
for the same four discs as in Fig. [I] As can be seen, /j, varies 
around a mean of 2 - 3. Hence, the values of pt roughly 
agree with the overall value of 2.6 and are comparable to 
the value of ~ 2 found in run 2.6-4. Moreover, p is well in 
the range where simulations without turbulent motions have 
found sub-Keplerian discs only. We therefore conclude that 
turbulent reconnection is not responsible for the build-up of 
Keplerian discs in our runs. 

Another way of reducing the magnetic braking efficiency 
was investigated by Hennebelle & Ciardi (20091 and Ciardi 



& Hennebelle (20101. These authors found that even for a 
small misalignment of the overall magnetic field and the ro- 
tation axis Keplerian discs can form. As we consider a tur- 
bulent flow, it is very likely that the magnetic field and the 
rotation axis are misaligned. In the right panel of Fig. [3] we 
plot the angle a between the disc angular momentum vec- 
tor and (B) in the spheres around the discs of run 2.6-4-A. 
The angle a is significantly larger than 0° which supports 
the picture of a reduced magnetic braking efficiency due to 
a misalignment of the magnetic field and the rotation axis. 

However, there is a third way to reduce the magnetic 
braking efficiency while simultaneously keeping the inwards 
angular momentum transport on a high level. Considering 
the top panel of Fig. [2] it can be seen that in the surround- 
ings of each disc there is a turbulent velocity field with no 
signs of a coherent rotation structure. Therefore no toroidal 
magnetic field component (w.r.t. the coordinate system of 
the disc) can be built up. But as the angular momentum 
is mainly extracted by toroidal Alfvenic waves, it is not 
surprising that the magnetic braking efficiency is strongly 
reduced in the environment of the disc despite a low mass- 
to-flux ratio (compare left panel of Fig.j3|. Despite the lack 
of a coherent rotation structure, locally the inwards angu- 
lar momentum transport can remain high due to local shear 
flows driving large angular momentum fluxes. We also note 
that the non-coherent flow cannot be efficiently slowed down 
by the magnetic field as it does in case of large-scale co- 
herent motions. This can be seen in our previous simula- 
tions (Sei fried et al.pOll I without initial turbulence. Here 
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Figure 4. Ratio of T gas to T mag for the discs found in run 2.6-4-A 
(same colours as in Fig. [3]) . and for the disc in run 2.6-4 without 
turbulence (violet line). 



the angular momentum is removed almost completely be- 
fore the gas hits the disc (see also bottom panel of Fig. [2j|. 
Hence, it is the shear flow generated by turbulent motions 
that leads to the build-up of Keplerian discs. 

To quantify this, we calculate the torques of the gas T gas 
and the magnetic field r mag in cylinders of variable radii and 
a total height of 40 AlQ around the CoM of each disc. The 
symmetry axis of the cylinders is determined by the angular 
momentum vector of the corresponding disc. In Fig. [4] we 
plot the ratio of r gas to T mag for the four discs found in run 
2.6-4-A. On average, r gas exceeds r mag by at least a factor 
of a few, i.e. angular momentum is transported inwards at a 
higher rate than it is extracted by the magnetic field. This 
is also observed for the discs in the other runs not shown 
here. We note that the gravitational torques not shown here 
are generally even smaller than r mag . The only exception 
occurs for disc 2 in run 2.6-4-A where a strong spiral arm 
has formed exerting a strong gravitational torque (see right 
panel of Fig. [2t. We also note that the strong fluctuations 
of r gas /r mag for disc 3 (green line) are due to strong per- 
turbations in its vicinity. Disc 4, for example, perpendicular 
intercepts the plane defined by disc 3 causing the drop at 
r ~ 90 AU. 

In general, however, there is a net transport of angu- 
lar momentum inwards resulting in the observed build-up 
of Keplerian discs. In contrast, for the non-turbulent run 
r gas is almost perfectly balanced by r mag . Furthermore, com- 
paring the absolute values of r between the turbulent and 
non-turbulent case shows that - although this comparison 
is somewhat crude - in general in the turbulent case r mag is 
reduced significantly whereas r gas remains comparable. 

This fits in the picture described above where the mag- 
netic braking efficiency is reduced due to the absence of a co- 
herent rotation structure and the angular momentum trans- 
port by the gas remains high due to local shears flows. Hence, 
this confirms our assumption that the turbulent disc envi- 
ronment is responsible for the build-up of Keplerian discs. 

We also made a rough estimate for the magnetic brak- 
ing time T mag = pdisk/penv ■ Z/vau in the discs themselves, 
calculating the mean density pdisk and the scale height Z of 
the discs as a function of radius. Assuming a density contrast 



2 This height is found to be reasonable by visually inspecting the 
discs in edge-on view. Furthermore, a variation of this value does 
not qualitatively change the overall picture. 



Figure 5. Scaling of the magnetic field for different times in run 
2.6-4-A as well as in run 2.6-4. 



Pdisk/Pcnv = 100 shows that for the discs T mag is on aver- 
age a factor of 5 - 10 above the characteristic timescale of 
inwards angular momentum transport I/disk/T" gas . However, 
we again emphasise that the magnetic braking efficiency has 
to be reduced in the region external to the discs in order to 
allow for the build-up of Keplerian discs. 



4 DISCUSSION 



Recently Santos-Lima et al. (20121 have reported the for- 



mation of Keplerian discs in turbulent, strongly magnetised 
low-mass cores. These authors, however, attribute this to 
the effect of turbulent reconnection ( |Lazarian fc Vis hniac 
1999) lowering the magnetic flux. As shown in the left panel 
of Fig. [3] we find the value of /x on small scales to be com- 
parable with the overall value, indicating that there is no 
magnetic flux loss. To further support this, we plot the scal- 
ing of the magnetic field with the density in Fig. [5] The ob- 
served scaling B oc p ' 5 is very similar to the non-turbulent 
case (run 2.6-4 in Seifried et al. 20111 and does not change 
significantly over time. Hence, we argue that in our case no 
significant turbulent reconnection occurs and that magnetic 
flux loss is not responsible for the formation of Keplerian 
discs as proposed by |Santos-Lima et al.| ( |2012[ ). Further com- 
parison to their work is not possible as the authors do not 
consider the scaling of the magnetic field. 

We note that in some cases /1 (Eq. [T| slightly increases 
with decreasing radius and eventually reaches values above 
10 at radii < 100 AU. Hence, it could be argued that tur- 
bulent reconnection is happening on these scales. Also from 
Fig.|5|some flux loss is apparent at densities above ~ 10 -11 g 
cm Typical scales associated with this density are of the 
order of 30 - 40 AU. However, on these scales (< 100 AU) 
the velocity structure seems to be relatively well ordered 
(Fig. [2J which makes turbulent reconnection unlikely to hap- 
pen. Hence, we rather attribute the flux loss to numerical 
diffusion. However, we again point out that, for enough an- 
gular momentum being available to form a Keplerian disc on 
scales of ~ 100 AU, already on larger scales (> 500 AU) the 
magnetic braking has to be reduced. On such large scales, 
however, no significant magnetic flux loss is observed (see 
left panel of Fig. [3}. Hence, whatever accounts for the flux 
loss on scales < 100 AU, it does not affect the mechanism 
postulated in this work. 

The large misalignment of the magnetic field of the 
spheres and the angular momentum of the discs of up to 90° 
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(see Fig. [3| also might weaken the effect of magnetic brak- 
ing as suggested by 
fc Hennebelle] ( [2010 



Hennebelle & Ciardi ( 2009 1 and Ciardi 



In their work, however, uniformly ro- 
tating spheres were considered which is clearly not the case 
here (see Fig. [2J . Therefore, it is hard to tell to what extent 
in our case the misalignment affects the magnetic braking 
efficiency. Our disc analysis shows that the misalignment is 
not the main reason for the decreased magnetic braking ef- 
ficiency. Here we see that the angular momentum is "gener- 
ated" locally by turbulent shear flows and that the classical 
picture of magnetic braking for coherently rotating structure 
breaks down leading to the build-up of Keplerian discs. 

Recently, |Krasnopolsky eTa l. (2011) have proposed that 
including the Hall effect can result in the formation of large- 
scale Keplerian discs. They claim, however, that a Hall co- 
efficient about one order of magnitude larger than expected 
under realistic conditions would be required. Also Ohmic 
dissipation fails to produce Keplerian discs larger than 
roughly 10 solar radii in the early evolutionary stages ( |Dapp| 
fc Basu|2011 Dapp et al.|2011 1, unless a strongly enhanced 
resistivity is used ( Krasnopolsky et al. 2010 1. Hence, it seems 
that all three non-ideal MHD effects (see |Mellon fe Li|2009} 
for the case of ambipolar diffusion) have a hard time in ac- 
counting for the formation of Keplerian discs. However, as 
we have shown, already for the ideal MHD limit Keplerian 
discs can form in strongly magnetised cores when turbu- 
lent motions are included. Therefore it seems that non-ideal 
MHD effects or turbulent reconnection are not necessarily 
required to avoid catastrophic magnetic braking. 

Earlier observations of individual Class sources found 
cases where no well denned, Keplerian disc were de- 
tected (e.g. Belloche et al.|20"02 \. These observations suggest 
that discs might form at later stages which would alleviate 
the magnetic braking problem. However, comparing more 
recent observations of samples of Class sources with de- 
tailed radiative transfer models suggests that the majority 
of them harbours well defined protostellar discs forming in 
the earliest stages ( j0rgensen et al.|2009||Enoch et al.|2011[ ). 
We note that recently Hennebelle et al. (20111 and 



Commergon et al. ( 2011 1 also performed simulations of high- 
mass turbulent cloud cores finding large-scale outflows which 
require the existence of protostellar discs supporting our ar- 
guments here. 



5 CONCLUSION 

We have performed collapse simulations of strongly magne- 
tised (/j, = 2.6, 5.2), 100 Mq cloud cores. A turbulent veloc- 
ity field was superposed on the uniform core rotation. We 
find that after an initial collapse phase of ~ 15 kyr discs with 
typical masses of ~ 0.1 Mq form. The discs are up to 100 
AU in size and have Keplerian rotation velocities, a result 
in strong contrast to previous simulations of strongly mag- 
netised cores containing no initial turbulence. We showed 
that our findings neither depend on the randomly chosen 
turbulence field nor on the adopted cooling routine. 

We suggest that the main reason for Keplerian disc for- 
mation is the turbulent surroundings of the disc. As there 
is no coherent rotation structure on scales of several 100 
AU, the generation of a toroidal magnetic field is suppressed 
therefore lowering the magnetic braking efficiency already 



before the gas hits the disc. At the same time the inwards 
angular momentum transport by the gas remains high due 
to local shear flows which results in a net inwards angular 
momentum transport. 

Our work strongly suggests that the "magnetic braking 
catastrophe" as reported in numerous papers is more or less 
a consequence of the highly idealised initial conditions ne- 
glecting turbulent motions. A turbulent velocity structure in 
the surroundings of the disc, as obtained with more realis- 
tic initial conditions, is enough to allow for the formation of 
Keplerian discs. Other effects like misaligned magnetic fields 
and angular momentum vectors, turbulent reconnection or 
non-ideal MHD effects seem not to be necessary. Turbulence 
alone provides a natural and at the same time very sim- 
ple mechanism to solve the "magnetic braking catastrophe" 
problem. 
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